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Abstract 

We study the matter produced in heavy ion collisions assuming that this matter is 
strongly interacting and employing AdS/CFT correspondence to investigate its dynamics. 
At late proper times r we show that Bjorken hydrodynamics solution, obtained recently 
by Janik and Peschanski using gauge-gravity duality [1], can be singled out by simply 
requiring that the metric tensor is a real and single-valued function of the coordinates 
everywhere in the bulk, without imposing any constraints on the curvature invariant. At 
early proper times we use similar strategy to show that the energy density e approaches 
a constant as r ^ 0. We therefore demonstrate that the strong coupling dynamics 
incorporates the isotropization transition in heavy ion collisions. By matching our early- 
time regime with the late-time one of Janik and Peschanski we estimate the isotropization 
time at RHIC to be approximately riso ~ 0.3 fm/c, in good agreement with results of 
hydrodynamic simulations. 
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1 Introduction 



The problem of thermalization of the system produced in heavy ion colhsions is crucial for our 
understanding of strong interactions. While in recent years there has been some significant 
progress in understanding the initial (pre-thermalization) stages of heavy ion collisions in the 
framework of the Color Glass Condensate (CGC) [2-19] (for a review of CGC see [20-22]), 
a complete theoretical understanding of the subsequent isotropization and thermalization of 
the produced medium continues to evade us. The success of hydrodynamics simulation in 
describing the data generated at RHIC is quite impressive: however, it requires a very early 
starting proper time for hydrodjTiamics, t<0.5 fm/c [23-30]. The existing theoretical models 
based on the weak-coupling dynamics so far have not been successful in reproducing such an 
early thermalization time [31-37]. Moreover, as was suggested by one of the authors [38-40], 
the perturbative weak-coupling dynamics, describable by Feynman diagrams, may in principle 
be unable to generate isotropization, and, therefore, thermalization in high energy heavy ion 
collisions. It was suggested in [38-40] that the only way isotropization and thermalization can 
be achieved in heavy ion collisions is due to strong-coupling effects. 

Another more phenomenological argument indicating that strong coupling dynamics could 
be important for heavy ion collisions at RHIC is the fact that hydrodynamic description of the 
collisions [23-30] has to have a very low shear viscosity [41]: the quark-gluon plasma produced 
at RHIC is almost a perfect fluid. Since the perturbative shear viscosity scales as r/ ~ 1 /(^^ [42] 
one concludes that low viscosity requires large coupling. Such conclusion was confirmed by 
explicit calculations employing gauge-gravity duality in [43-47]. 

Unfortunately, the strong coupling regime of QCD is not under theoretical control. However, 
for some gauge theories, like A/" = 4 super Yang-Mills theory, the strong coupling regime can be 
studied using AdS/CFT correspondence [48-51]. To understand the strong coupling dynamics 
of the matter produced in heavy ion collisions one has to model QCD with A/" = 4 super 
Yang- Mills theory and use AdS / CFT correspondence to study the strong-coupling dynamics of 
such theory. In recent years there have been many applications of AdS/CFT correspondence 
to heavy ion collisions [1,52-62]. 

The goal of this paper is to attempt to understand the onset of isotropization and ther- 
malization in the AdS/CFT framework. We will consider an idealized Bjorken picture of the 
collision [63], in which nuclei have an infinite extent in the transverse plane and the distribution 
of the produced matter is rapidity-independent: the energy and pressure components of the 
energy-momentum tensor in the local rest frame would then depend only on the proper time 
r. We would like to emphasize that isotropization and thermalization of the produced medium 
do not necessarily happen at the same time: isotropization is necessary for thermalization, 
but not vice versa. As was argued in [14], in the CGC framework the energy density e of the 
matter produced in heavy ion collisions scales as a power of a logarithm of the proper time r, 
as shown in Eq. (1551) below. Therefore, the produced matter is anisotropic: its longitudinal 
pressure component is in fact negative. Isotropization for such matter would imply generation 
of positive longitudinal pressure, comparable to the transverse pressure components. Note that 
for isotropization to take place it does not appear necessary for the system to have a thermal 
distribution: hence thermalization may take place after isotropization, and not necessarily at 
the same time. Below we will study the onset of isotropization, since it is much easier to infer 
from the components of the energy-momentum tensor. 
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In a pioneering paper [1] Janik and Peschanski showed that for the strongly-coupled medium 
produced in heavy ion collisions and described by gauge-gravity duality the late time asymp- 
totics can only lead to Bjorken hydrodynamics [63] with the energy density scaling with proper 
time as e ~ r~^/^. Since then a number of further investigations have been carried out [53,54,64], 
mostly concentrating on the subleading at late times viscous corrections to the ideal hydrody- 
namic behavior of [1]. In this work we will take a different route and concentrate on the early 
time behavior of the strongly coupled medium produced in heavy ion collisions. The hope is 
to verify if such medium is anisotropic at early proper times: if that proves to be the case 
then the onset of isotropization would necessarily be a part of strongly-coupled medium's dy- 
namics at intermediate proper times. The idea of studying the matter produced in nuclear 
collisions at small proper times has been originally suggested in the CGC framework in [14,65]. 
Here we will apply the same concept to the strong-coupling dynamics employing the AdS / CFT 
corresp ondence . 

The paper is structured as follows. We begin in Sect. [2] by setting up the formalism of 
holographic renormalization in the context of heavy ion collisions with Bjorken geometry. We 
then proceed in Sect. [3] by solving Einstein equations ([3]) in the AdS space assuming that the 
four-dimensional energy density scales as a power of proper time, as given by Eq. (fT7|) . In Sect. 
13. II we will demonstrate the first few steps of the perturbative solution expanding the metric in 
powers of the Fefferman- Graham coordinate z [66]. In Sect 13.21 we proceed by re-analyzing the 
late proper time dynamics: we show that in order to pick the correct Bjorken hydrodynamic 
behavior of the solution of Einstein equations ([3]) one only has to impose the condition that the 
resulting metric is a real single- valued function of the coordinates everywhere in the bulk. Hence 
we demonstrate that the condition of non- singularity of the curvature invariant employed in [1] 
is not needed to constrain the dynamics in four dimensions from AdS/CFT. We then apply 
the same approach to study the early time dynamics in Sect. 13. 3t we solve Einstein equations 
([3]) at small-r and, after imposing the condition for the metric to be real and single- valued, we 
derive that energy density can only go to a constant value as r ^ (see Eq. flSTl) ). 

This is our main result: as the system having constant energy density in heavy ion collisions 
also has negative longitudinal pressure, we conclude that at very early times the strongly inter- 
acting system described by gauge-gravity duality is very anisotropic. At the same time we know 
that at late proper times it becomes isotropic evolving according to Bjorken hydrodynamics. 
We thus demonstrate that the isotropization transition does take place at intermediate times 
for the evolution of this strongly-coupled system! Moreover, while indeed the early-time energy 
density in CGC [Eq. (15^ ] and in the strongly-coupled medium [Eq. (IHT]) ] behave somewhat 
differently, the two initial-time behaviors are not entirely unlike each other: in both cases the 
initial longitudinal pressure is negative, equal to the negative of the energy density, while the 
transverse pressures are equal to the energy density [see Eq. fl52l) ]. We therefore conclude that 
the behavior of the strongly-coupled system studied here may be relevant for real-life heavy ion 
collisions, where the initial stages of the collisions are successfully described by the CGC. 

We verify our result of constant energy density at early times in Sect. IH where, inspired by 
CGC [14], we look for possible logarithmic scaling of the energy density at early time, but find 
none. The energy density remains constant at small r, without any logarithmic divergence. 

In Sect. |5] we match the early time dynamics considered in this work onto the late-time 
Bjorken hydrodynamic solution found in [1] to estimate the isotropization time (see Fig. [T]). 
The result is given by Eq. (15T]) . If we apply this formula to central Au -|- Au collisions at RHIC 
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at ^/s = 200 GeV we obtain the isotropization time of Tiso ~ 0.3 fm/c, in good agreement with 
results of hydrodynamic simulations. 

We conclude in Sect. Elby restating our main results and by summarizing our knowledge of 
energy density as a function of proper time for the strongly-coupled system at hand in Fig. [2j 

2 Holographic renormalization: general setup 

Let us consider a head-on collision of two very large identical nuclei. For simplicity we assume 
that the distribution of matter produced in such collisions is independent of (space-time) ra- 
pidity r]. Since the nuclei are identical the matter distribution should also respect x+ x_ 
(?7 ^ —f]) symmetry. Finally, as the nuclei are very large, we assume that they have infinite 
extent in the transverse directions. Therefore the energy-momentum tensor of the produced 
matter is independent of the transverse coordinates x = (xi,X2). This is indeed the same 
approximation as employed by Bjorken in [63]. 

Following Janik and Peschanski [1] we write the metric in AdSs space for such a system 
using Fefferman-Graham coordinates [66] as 

ds'^ = 4 dr^ + B(t, z) drf + C(r, z) dx\ + dz'^] (1) 

z^ 

where r = a/2x+ X- is the proper time in four dimensions, rj = (1/2) ln(x+/x_) is the space- 
time rapidity, x± = (xq ±0:3) /a/2, dx^ = dxf + dxl and z is the coordinate in the 5th dimension. 
Just like in [1] we write 

A{t,z) = e'^("'"), B{t,z) = e''^^''\ C{t,z) = e^^^''\ (2) 

where a(r, z), b{T, z) and c(r, z) are some unknown functions, to be determined from the solution 
of Einstein equations with a negative cosmological constant A = —6 in AdSs space 

Rfiu -^QfiuR-Ggfiu = 0. (3) 
Here i?^,^ is the Ricci tensor and R is the scalar curvature. Contracting Eq. ([3]) with g'^'^ yields 

R = -20 (4) 
which, when inserted back into Eq. ([3]) gives 

Rf^u + 4^g^^ = 0. (5) 

Finding an exact general solution of Eq. with the metric ansatz ([T]) appears to be a 
daunting task. Instead, we begin by exploring the solution using the technique of holographic 
renormalization [67]. The general prescription for a metric of the form 

ds^ = ^ [g^^{x,z)dx^dx'' + dz^] (6) 

is to expand it near the ^ = boundary in a series [67] 

-qA^.z) = ~g^^J{x)+z'~gl^Xx)+z'~g;^Xx) + ... (7) 
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and solve Einstein equations ([5]) order-by-order in z. 

In our case we know the metric on the 2 = boundary: it is simply the fiat Minkowski 
metric, such that gliJ{x) = rj^y. Then one can show that gliJ{x) = [67]. Finally, the fourth 
order term in Eq. ([7]) is related to the expectation value of the energy-momentum tensor T^j, 
in our four dimensions [53, 67] 

]\t2 

{T,u) = ^ ^it^(x) (8) 

with Nc the number of colors. As can be easily shown [38] , the most general energy- momentum 
tensor in the boost-invariant geometry of the collision of two very large nuclei can be written 
(in the co-moving frame) as 

/e(r) \ 
p{t) 
p{t) 
V p3{r)J 



(9) 



in the t, xi, X2, X3 coordinate system. Here e(r) is the energy density, p(r) is the transverse pres- 
sure component, and P3(r) is the longitudinal pressure component. (The latter two quantities 
do not have to be equal in general: they are equal only in the case of ideal hydrodynamics.) 

In [1] it was elegantly shown that requiring non-negativity of the energy density e(r) in all 
frames leads to the following conditions on its dependence on the proper time r 

e'(r) < 0, re'(r) > -4e(r). (10) 

We will use this conditions later to restrict energy density. 

As will be manifest shortly, Einstein equations lead to the following conditions for the 
energy-momentum tensor. First of all they require energy-momentum conservation in four 
dimensions 

d, = (11) 

which, when applied to Eq. ([9]), gives 

^ = -i±^. (12) 

dr T 

(Here and below we will sometimes suppress the arguments of e, p and p^.) Secondly they 
require the energy-momentum tensor to be traceless 

(t;) = 0, (13) 

which, for Eq. leads to 

e = 2p + p3. (14) 

Finally, to apply the expansion of Eq. with gj^J (x) = ?7^i, and g^S (x) = to the metric 
in Eq. ([1]) we write [1] 

00 00 00 

a(r,z) = ^a„(r)/+2^ b{T, z) = J^bnir) z'+'\ c(r, z) = ^ c„(r) 2^+^". (15) 

n=0 n=0 n=0 
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Comparing these to Eqs. ([8]) and ([9]) yields 



Ar2 jv^ iV2 

ao(r), P3 = tH? boir), p = Co(r). (16) 



Now we are ready to implement the holographic renormalization using a particular ansatz 
for the energy density. 

3 Power-law scaling of energy density with proper time 

3.1 First steps: iterative solution 

Let us first assume that the energy density scales as 

e(r) ~ r^. (17) 



This is the simplest assumption. We know that it works at late times r ^ 1 (in units of yGp^ 
with Gtv the Newton's constant in AdSs) [1]. In a little while we will try to apply it to the 
early-time dynamics (r <^ 1) and if it would not work we will be forced to look for a more 
sophisticated one. At the moment we are not going to make any assumptions or approximations 
specific to early or late times. 

We will begin with the scaling of Eq. (fT7|) and try to build a solution of Einstein equations 
in the bulk using holographic renormalization. From the dynamics in four dimensions we know 
that if e ~ then p ^ and ~ r^. Using Eq. ( |T6l) we write 



ao(r) = aoT^, boir) = boT^, co(r) = cqt^ (18) 

with ao, bo and Cq some constants. Substituting Eq. ffTSl) into Eq. f|T5|) . and using the expansion 
from the latter in Eq. ([5]) we solve Einstein equations to the lowest non-trivial order in z. This 
gives 

ao(r) = aoT^ (19a) 
boir) = ao(A + l)r^ (19b) 

Co(r) = -ao — —T . (19c) 

Indeed Eqs. (1191) could have been obtained from Eqs. ( |T2l) and ( |T4|) directly from the four- 
dimensional considerations of dynamics without trace anomaly. Here we obtain them from the 
AdS/CFT correspondence. 

Substituting Eqs. (fT9!) back into Eq. (fT5|) . and again solving Einstein equations ([5]) at the 
lowest non-trivial order in z (which is now different from the previous step) yields 

a,(r) =a„^^plT^--' (20a) 
,,M.«„^(^±|(^,- (20b) 

ci(r) = -ao \^ ' r^-\ (20c) 
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For what we intend to do next it is instructive to iterate the procedure several more times. 
Here we will show only the result of the next iteration: 



a2(r) = --^[aor^-'(4A2-A^)+a2r2A8(8 + 8A + 3A2)] (21a) 

b^(r) = [aor^"'(-12A2 + 4A-V3A'^-A5) + ag r^^ 8(8 + 8 A + SA^)] (21b) 

C2(r) = --^ k r^"' (8 A^ - 4 A=^ - 2A^ + A^) + 16 (8 + 8 A + 3 A^)] . (21c) 

768 

As follows from the conditions in Eq. (ITU]) [1], positivity of energy density e combined with 
the power-law ansatz of Eq. fll7p requires that 

-4 < A < 0. (22) 

Let us assume that A > —4. This is a safe assumption: at late times the smallest possible 
value of A is A = —4/3 [63]. At late times we therefore have —4 < A < 0. 

At early times A > — 1, the total energy of the matter produced in a collision approximately 
scales as ~ e r for small r. For A < — 1 the total energy E would be infinite in the r — ^ limit, 
which is impossible. Hence the physically relevant values of A at early times are constrained 
to -1 < A < 0. 



3.2 Re-deriving Bjorken hydrodynamics 

Let us first re- visit the dynamics at late times, which has already been studied in [1]. If A > —4, 
then the r^^ term dominates over r^~^ term in Eq. (12T1) . such that Eqs. ( fT9l) . (l20l) and (l2Til 
combined with Eq. f|T5l) give us the following type of a series for a(r, z) in the large-r limit 

a(r,z) = #z^r^ + #'z^r2^ + ... (23) 

with # and some A-dependent constants. (The series for 6(r, z) and c(r, z) are similar.) 
One can see that a(r, 2) becomes a function of a single scaling variable v oc zr^^^, as was 
observed and utilized in [1]. 

Following Janik and Peschanski [1] we define the scaling variable as 

V = (-a^y/^zT^/^ (24) 

Note that combining Eqs. 0161] and (ITUh ) we get 

Z 7r^ 

such that positivity of the energy density requires that < and the power of — ao in Eq. (12^ 
is well-defined. Assuming that a(r, z) = a{v), b{T, z) = b{v), and c{t,z) = c{v) the authors 
of [1] solved the Einstein equations ([5]) for the metric of Eqs. ([1]) and ([2]), keeping v constant 
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and taking r oo limit. This led to the following solution obtained in [1] 

a{v) = l{l-^) HI + Dv')+^-(^l + ^yn{l-D v') (26a) 

h{v) = ^ (l - ^) + D v') + ^ (l + ^) - D v') (26b) 
1 , A + 2\ , ^ 4, 1 A A + 2\ 



where 



'^''^ = 2 + ^''^^ + ^''" + 2 (1-^^) Ml-^^'), (26c) 



To determine the allowed values of the power A let us note that ln(l — Dv'^) has a branch 
cut along the real axis for 1 — Dv^ Generally speaking this means that A(r, 2;), B{t,z) 
and C(r, z) will be complex- valued for 1 — D v"^ < 0. Since the t —>■ 00 limit taken in arriving at 
Eq. ( 1261) was done keeping v fixed there is no constraint prohibiting v > D^^^^ and preventing 
the metric from becoming complex. The only way for the metric to be real and single-valued 
for all V is if the coefficients in front of ln(l — Dv^) in Eq. (l26l) are integers. We therefore 
require that 

^ 1 + ^ ) = n (28a) 



2 V D 

^ 1 + ^-r^ ] = m (28b) 



2 V D 

1 f 1 - ^) = , (28c) 

2 \ 2D J ^ ' 

with n, m, I some integer numbers. After some straightforward algebra one derives 

n + m = 2(1 - /). (29) 

Therefore n -|- m is an even number. Then n — m is also an even number. We then write using 
Eq. (IMd) 

^ n-m = 2k (30) 



2D 

with k some other integer. Solving 

A = -4:Dk (31) 

for A with the help of Eq. fl27p yields 

4 k 



^ ~~ 3 2 A;2 _ 1 L 



2 A; ± V3 - 2 A;2 



(32) 
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As A < and D > we conclude from Eq. (13T]) that k > 0. For the solutions in Eq. (!32l) to 
be real, the only allowed non-negative values of k are k = 0,1. 

li k = then A = and D = 1/ ^/3. However, then Eq. fl28b ) will not be satisfied as n will 
not be integer. 

li k = 1 the allowed values of A are A = —4 and A = — |. However, as we concluded from, 
say, Eq. fl2T|) . for A = —4 the solution of Eq. fl26l) is no longer dominant, as terms subleading 
for larger values of A become comparable to it. On top of that A = —4 is not a viable physical 
value of the power. 

We are then left only with A = — |, which, when used in Eq. ( 125|) . gives us 

iV^ 1 1 

= -^«o^ (33) 

characteristic of ideal Bjorken hydrodynamics [63]. We have thus obtained the perfect fluid 
behavior for the matter produced in heavy ion collisions without imposing a constraint of the 
absence of singularities in the curvature invariant TZ = Rn^pcr R^'^^^ , like it was done in the 
original derivation of [1]. The only necessary requirement was that the metric tensor is a real 
and single-valued function of v for all values of the scaling variable v. 



3.3 Solution at early times: scaling ansatz 

We will now proceed towards the main goal of this work, which is studying the early-time 
asymptotics, r ^ 1. For small r the dominant part of the holographic renormalization series 
studied in Sect. 13.11 is 

a(r, z) = ifz^T^ + #' ^ .8 ^ _ _ _ 

= ^'r^ (# + #'^ + #"^ + ---) (34) 

with 7^, and 7^" some (different) A-dependent constants. (The series for 6(r, z) and for 
c(r, z) have the same structure.) Unfortunately, no single scaling variable can describe the r 
and z dependence now. Nevertheless, as one can see from Eq. (1341) . the series is in z/t. Hence, 
defining a new scaling variable 

u=- (35) 
r 

we rewrite Eq. as 

a{T,u) = T^+^u'' + + + ...). (36) 

Therefore, a general ansatz for a(r, m) at early times is 

a(r,u) = T^+^a{u) (37) 
where a{u) is some unknown function. Similarly we write 

b{T,u) = r^+'/5(u), c{t,u) = r^+S(«)- (38) 
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The initial conditions for a{u), P{u) and 7(m) can be determined from Eqs. (fT9!) yielding 

A + 2 

a{u) = aou'^, P{u) = ao (A + 1) u^, 7(u) = — — - — u"^, as u ^ 0. (39) 

Our goal now is to use the ansatz of Eqs. (137|) and (l38l) and solve Einstein equations ([5]) in 
the r — s> limit keeping u fixed. In this limit we write 

A{t,u) = e"("'") = e"^^'°(") = l + r^+"a(M) +o(r2^+S) (40a) 
B{t,u) = l + T^+^f3{u) +o(r2^+«) (40b) 

C(r,n) = 1 + +o(r'^+') (40c) 

as, with the ansatz of Eqs. (|37j) and (|38l) . we do not have any control over o(r^^^^) terms in 
A, B and C. 

Rewriting the metric ([T]) in terms of r, u, r/, and x_ = {xi,X2) we plug it into Eq. ([5]). At 
the lowest order in r we obtain the following equations, corresponding to the uu, u r and xiXi 
(or, equivalently, X2X2) components of the Einstein equations: 

a' 3' 1' a" B" 

— + ^ + ^- — - ^-7" = 0, (41) 
a' - ^ /5' - (2 + A) y - ^ n a" = 0, (42) 



+ 4A + -A2 7 + ^ + ^+ h^-^^-An 7' + -(n^-l)7" = 0. (43) 
2 J 2u 2u \2u 2 J 2 

Here a' = da{u)/du, a" = a{u) / dv^ , and the same applies to derivatives of I3{u) and 7(u). 
After an integration using the initial conditions ( l39l) . Eq. (jHj) gives 

a + /5 + 27 = 0. (44) 

Using Eq. (jH]) we eliminate 7 from Eq. ( I42l) to obtain 

/5 = (5 + A)a-ua'. (45) 

Substituting Eqs. (iS]) and dHj) into Eq. (iSl) yields 

{u^ - 1) a'" + n [7 + A - (11 + 3 A) n^] a" 
+ [-3 (5 + A) + (51 + 25 A + 3 A^) u^] a' - (4 + A)^ (6 + A) n a = 0. (46) 

The solution of Eq. fHB]) satisfying the initial condition (jSHD can be easily found by searching 
for it in the form of an infinite series in u starting at the order u^. From Eq. fl5B]l we also 
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infer that the coefficient in front of the term in the series is zero. The solution of Eq. (l46l) 
satisfying such conditions is 



a[u] 



F 



A A 



(47) 



with F the hypergeometric function. Eqs. (HSj) and (1441) allow us to find (3{u) and 7(m). In the 
end, using Eqs. (HOj) . we obtain the following solution for the components of the metric ([1]) 



A{t,u) 



B(t.u) 



1 + aQT 
1 



4+A 



u'^F -1 



A 

2"' 



A 



4+A 4 



(A + 1)F -1 



A A 



(48a) 



2 ' 



C{t,u) 



4+A 4 

ao r u 



A (A + 2) 
6 

A + 2 
12 

+ Am^f 



A A 



M^F 1 , ;4;n 



A 



(48b) 



A 



A 


A 






;4; 





(48c) 



As is well-known, hypergeometric functions, such as we have in Eq. (1481) . have a branch-cut 
discontinuity along the real axis for m > 1. Therefore, for the metric to be a real and single- 
valued function of u at u > 1 we need the hypergeometric series to terminate. The reasoning 
here is the same as what we used in re-deriving Bjorken hydrodynamics in Sect. 13. 2[ As was 



pointed out at the end of Sect. 13. ![ the allowed values of the power A at early times are 
— 1 < A < 0. Within this range it is only possible for the hypergeometric series to terminate if 
A = 0. Therefore, A = is the only allowed solution! 
Plugging in A = into Eq. (148|) gives 



A{t, u) = 1 + a^T^ 
B{r,u) 
C{t,u) 

Using A = in Eqs. ^ and ^ yields 



1 I 4 4 

1 4 4 



«0, P3 



Nl 
2 7r2 



ao, p 



(49a) 
(49b) 
(49c) 



(50) 



Note that at early times ao is different from the one used at late times in Sect. 13.21 Still 
ao < and we conclude from Eq. fl50|) that, for the strongly-coupled medium described by 
gauge-gravity duality in the geometry of heavy ion collisions 



We also obtain that, as r 



e(r) - 




constant 



as 



0. 



p{t) = e(r). 



(51) 



(52) 
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The latter relations one can also obtain from Eqs. (fT2i) and (fT4|) by using Eq. (!5T!) in them. 

Note that for the description of heavy ion collisions in the framework of the Color Glass 
Condensate (CGC), in which the produced medium has a small coupling constant, it was 
suggested recently by Lappi [14] (see also [68]) that the leading behavior of energy density as 
r ^ is given by 

ecGc(i") oc In^ as r ^ 0. (53) 

Using Eq. (l53l) in Eqs. (fT2l) and iHM and keeping leading logarithms only one arrives at Eq. (l52l) . 
Hence, while Eq. (1531) is indeed different from our Eq. fl5T|) . the energy density and pressures 
at early proper times behave somewhat similar for the cases of CGC and the strongly coupled 
medium considered above. In both cases the initial longitudinal pressure is negative, while the 
transverse pressure is positive and equal to the energy density. 



4 Logarithmic scaling of energy density 
with proper time 

In the previous Section we looked for the energy density at early proper times using the power- 
law ansatz (fTTIl and found that the power of r has to be zero. However, in view of the CGC result 
(1531) [14] it may be that the absence of a power-law scaling indicates some residual logarithmic 
scaling of energy density with proper time. Here we will verify whether such logarithmic scaling 
does indeed take place. 

It is tempting to begin by assuming that the energy density scales logarithmically with r 
such that 

e(r) oc In^ (-] (54) 



with some number 6 > 0. However, if we try solving Einstein equations ([5]) with the ansatz ([5i 
at z = keeping the leading logarithmic [ln(l/r)] terms at small r we will simply reproduce the 
results of the previous Section [with the second terms in each of the equations fHOj) multiplied by 
In"^ (^)] without obtaining any constraints on 6. The reason for that is quite straightforward: 
if one differentiates a power of ln(l/r) with respect to r one power of the logarithm would 
be lost and such contribution would not be leading logarithmic anymore and would have to 
be discarded. Therefore, in the leading logarithmic approximation no derivative acts on the 
power of the logarithm, making In'' (i) term just an overall factor canceling out in the Einstein 
equations and not giving any constraint on S. 

To avoid this problem one has to include the subleading logarithmic correction. We start 
by writing 

ao(r) = ao In^ f - ) + ai In^'^ (-)+••• (55) 



such that, using Eq. (jT6l) . 

e(r) 



2 7r2 



r / \ r 



ao In^ ( H + «i ( ^ 



(56) 
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Here Oi is some undetermined constant. The ellipsis denote terms of the order o [in'' ^ (~)]- 
Using Eq. dSSD in Eqs. and ([HD yields 



6o(r) = ao In^ (^^j + {a, - 5) In^"^ (^j + " " " 

Co(r) = -ao In^ (i) + " «i) 1^'"' (7) + " " " " ^^^^^ 

We performed a perturbative expansion of the solution of Einstein equations, similar to the 
one carried out in Sect. 13.11 It indicated that the scaling variable is still u defined in Eq. ( !35l) . 
Therefore we try looking for the solution of Eq. ([5]) in the following form 

A{t, u) = 1 + (^^^ + ai{u) ln^~' (7) + ' ' ' ^^^^^ 

B{t, u) = 1 + \Yi^ + I3i{u) In^-^ (7) + " " " ^^^^^ 

C(r, u) = l-aoT^ In^ + 71 (m) In^"^ (7) + ' ' ' ^^^"^^ 

with the unknown functions ai{u), (3i{u) and 71 (u) satisfying the initial conditions which follow 
from Eqs. ([55]), ([57D 

ai(M) = aiM^, /3i('u) = (ai — ao^)w^, 71 (m) = ^-^ cti^ m^, as m — >• 0. 

(59) 

The ansatz of Eq. fl58|) trivially satisfies the Einstein equations at the leading logarithmic 
order for the reasons discussed above. Non-trivial information comes in at the sub leading 
logarithmic order. At the order In'^"^ (i) the nn, xiXi and uxi components of the Einstein 
equations (jS]) read 

«'i I3\ 7i' «'/ „ . ^ 

/ n 2 u n 2 2 

16M7i + a;+/3i + (5-7M2)7; + ^i(M2_^)^// _ ^g^^ 

Aa^Su^ -2a[ + ?>(3[ + A-i[ + ua'{ = 0. (62) 
Just like in Sect. 13.31 we integrate Eq. (!60l) to obtain 

«!+/?! + 271 = 0. (63) 
Using Eq. (163|) in Eq. (16T|) to eliminate a[ + yields 

16m7i + (3-7m2)7; + m(m2-1)7;' = 0. (64) 
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The only solution of Eq. (IMI) satisfying the initial conditions (159|) is 



7i(n) = [^-a,)u\ (65) 



2 

Using the solution (!65ll in Eqs. (l62l) and (j63ll yields 

M < - 5 + 8 ai = 0. (66) 
The solution of Eq. fl66l) satisfying the initial conditions fl59|) is 

«!(«) = Oi M^. (67) 

Finally, plugging Eqs. (!65|) and (167|) into Eq. (163!) yields 

/5i(m) = {ai — ao6)u^. (68) 

To make sure that Eqs. (1671) . (l68l) and (l65l) give us the solution of Eq. (I5l) we need to verify that 
the remaining two equations corresponding to 7777 and rr components of Eq. ([5]) are satisfied. 
These read 

2aod - Auai + 20u pi + 8u-fi + {1 + u'^) a[ + 4:{1 - 2u^) p[ + 2 {1 - u^) 

-u (1 - M^) p'l = (69) 

and 

-6 ao + 4 n ai - 20 n /?i - 24 n 7i + 2 (n^ - 2) - (n^ - 1) (-/^^ - 2 7^ + n a'/) = 0. 

(70) 

Substituting Eqs. ( 1671) . ( !68l) and ( l65l) into Eqs. (!69l) and (1701) we see that the latter equations 
can be satisfied for Cq 7^ only if 

6=0. (71) 
This constrains the power of the logarithm in Eq. ( l54l) . Putting 5 = in Eq. ( 156|) we obtain 

^' as r^O. (72) 



e r 



2 7r2 



ao + 



ln(l/r) • • 

We therefore confirm the result of the previous Section that the energy density approaches a 
constant as r — 0. Eq. (172|) demonstrates that the approach to the constant may be logarithmic. 

Again our derivation never employed the requirement of no singularities of the curvature 
invariant TZ = Rfj^upa R^^'^" of [1] . At the same time it can be explicitly checked that it does not 
violate such a requirement. The curvature invariant for the solution in Eq. (H^ (which is also 
the leading small-r solution for the logarithmic scaling considered in this Section, as follows 
from the above) is 

7^ = 40 + o(r«). (73) 

Here we keep u fixed and take r to be small. Eq. (173|) shows that the non-trivial part of the 
solution (149|) does not contribute to the curvature invariant and, therefore, does not introduce 
singularities in it. In fact one can show in general that the first perturbative correction to the 
metric of the fiat AdS space satisfying Eq. ^ does not contribute to the curvature invariant. 
It appears that it would have been impossible to constrain the early times scaling of the energy 
density using non-singularity of the curvature invariant only. Luckily, we were able to pinpoint 
the right metric without imposing any additional conditions. 
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5 Estimate of isotropization time 



We have shown above that the strongly-coupled system produced in heavy ion collisions with 
the dynamics described by the AdS/CFT correspondence starts off at early times with a finite 
energy density and negative longitudinal pressure component of the energy-momentum tensor. 
Hence, the system at these early times is very asymmetric. On the other hand, as was demon- 
strated by Janik and Peschanski [1], the same strongly-coupled system at asymptotically late 
times behaves like the ideal Bjorken hydrodynamics, with the longitudinal and transverse pres- 
sure components equal and positive. Therefore, the dynamics of this strongly-coupled system 
incorporates the isotropization transition at some intermediate time, and has the potential of 
solving the thermalization problem in heavy ion collisions. Indeed, since the initial condition in 
real heavy ion collisions described by QCD are likely to be perturbative describable by CGC, 
to verify the relevance of the strongly-coupled model considered above to the real-life heavy 
ion collisions it is important to check whether the early times dynamics obtained here agrees 
with perturbative expectations. The very early time dynamics was studied in the CGC frame- 
work in [14] with the conclusion that the energy density scales as shown in Eq. ( l53l) . Such 
behavior also leads to negative initial longitudinal pressure, similar to what we observed for 
the strongly-coupled theory. Therefore it is plausible that the strongly-coupled isotropization 
scenario is relevant for the real- life heavy ion collisions described by QCD. 

With that in mind let us estimate the onset of isotropic Bjorken hydrodynamics in the gauge- 
gravity duality framework considered here. First of all, at early times the general solution PHI) 
has a branch-cut discontinuity at m > 1. In the 2, r-space the branch cut is at 

z > zl'^'^y = T. (74) 

On the other hand, the solution found by Janik and Peschanski [1] and shown here in Eq. fl2Bl) 
has a branch cut at 

V > D--i (75) 
which, with the help of Eq. (l24l) with A = —4/3 and corresponding D = 1/3 translates into 

. > ^ ' . (76) 

The locations of the starting points of the branch cuts z'^J^^^^ and z^^^^ given by Eqs. fl74l) 
and (1751) are shown in the 2, r-plane in Fig. [H The two curves in Fig. [T] intercept at some time 
which we label Tiso- 

It is important to point out that, as can be seen, say, from Eq. f l21l) . the full solution with 
the power-law ansatz of Eq. (fT7|) consists of a sum (for a, h and c) of the solution in Eq. (!26|) [1] 
and the solution in Eq. (148|1 found above. While different solutions dominate at different times, 
both exist at all times. We now observe that in order to pinpoint the correct solution using the 
power-law ansatz of Eq. (IT7|) as was done above one starts outside the branch cut, i.e., aX z < Zc 
for both early and late times, and imposes the no-branch-cut condition as z increases past Zc- 
For the times before the intersection of the two curves in Fig. [H r < Tiso, the first branch cut 
which has to be eliminated is the one starting at z'^^'^^^ . That fixes the power A = making 
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Figure 1: The locations of branching points z^^"^^^ and z^^^'^ for the early- and late-time metrics 
in the r-plane. The intersection of the two curves gives our estimate of the isotropization 
time Tiso (see text). 



the energy density a constant, as described above and shown in Eq. fl5T|) . At late times the 
first branch cut to be eliminated is at fixing A = —4/3 as described in Section [3l2l Hence 
the behavior of energy density, or, equivalently, the dominance of either one of the solutions 
in Eqs. f l2^ and fHSjl . is determined by which one of the two branch cuts starts at smaller 
values of z. At z^^"^^^ < z]?'^'^ the early-time solution dominates leading to the energy density 
scaling of Eq. (ISTl) . At z^^^^^ > z'^^'^ the late-time solution dominates leading to Bjorken's 
energy density scaling of Eq. (!33|) . While, strictly speaking both solutions are valid only at 
either asymptotically early or late times far away from the intersection point shown in Fig. [H 
we can estimate the time of the transition between the two regimes by the time of the intercept 
Tiso, defined by the condition 

yearly _ ^late_ ^77^ 

Using Eqs. ([74]) and (JTEl) in Eq. ((771) yields 

-) ' • (78) 

This is our estimate of the isotropization time, the time necessary for hydrodynamics to work 
and required for thermalization. To determine Oq we write Bjorken energy density scaling as 

<r) = ^. (79) 
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Comparing Eqs. (fTOl) and (!33l) we see that 



ao = 




(80) 



which, when substituted into Eq. (1751) gives 



( 



3 iVni 



(81) 



ISO 



Note that for a self-consistent conformally-invariant strongly interacting theory the above result 
is natural: from the scaling of Eq. (179|) we see that such theory would be characterized only 
by one dimensionful parameter — eg. Using this parameter we can construct the time scale, 
which is Cg'^^^. Indeed to find the coefficient in front of Bq^^^ in Eq. (IHT]) one needs to perform 
a more complete calculation, as was done above. 

To evaluate the isotropization time in Eq. (15T]) for central Au + Au collisions with ^/s = 
200 GeV at RHIC we use the fact that hydrodynamic simulations of [23-27], which successfully 
describe a variety of RHIC observables, yield the averaged over all impact parameters energy 
density of e = 15 GeV/fm'^ [27,69] at the proper time r = 0.6 fm/c. From these data, using 
Eq. ( !79l) we obtain eo ~ 38 fm~^/^. Substituting this number into Eq. ( ISTl) with Nc = 3 yields 



which is quite close to the initialization time for the hydrodynamic simulations required to 
describe RHIC data [23-30,69]! 

Indeed the isotropization time estimate of Eq. ( IHTi) depends on the final state observable eo 
instead of some initial-state observable. Therefore one should indeed question the predictive 
power of such estimate. We leave it for the future work to improve on this result and to tie it 
to some initial state observables. 

6 Conclusions 

We conclude by restating our main results. We have re-derived the result of Janik and Peschan- 
ski [1] stating that at late proper times the system produced in heavy ion collisions in the 
strongly-coupled regime described by gauge-gravity duality exhibits energy density scaling with 
time characteristic of Bjorken hydrodynamics, as shown in Eq. fl33l) . In such regime all three 
pressure components of the energy momentum tensor are equal to each other, as is required for 
a thermalized medium. To obtain this result we use a much weaker condition than that used 
in [1]: we require the metric to be a real and single- valued function of coordinates in the bulk. 

We then explore the dynamics of the same strongly-coupled system at early times and show 
that the energy density approaches a constant as r — > 0, as shown in Eq. ( l5Ti) . The longitudinal 
pressure of this early system is therefore negative [see Eq. (fT2|) ]. We have therefore demonstrated 
that this strongly-coupled system is initially anisotropic, but at late times it becomes isotropic 
and evolves according to Bjorken hydrodynamics. Therefore, a complete solution for the energy 
density of this system at all times would contain the isotropization transition! The similarity 



ISO 



0.29 fm/c 



(82) 
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Figure 2: The scaling of energy density e with proper time r in heavy ion colhsions for a strongly- 
coupled system as described by AdS/CFT correspondence. At late time Bjorken hydrodynamics 
is recovered, as was shown in [1]. At early times the energy density goes to a constant, as shown 
in this work. Our estimate of the isotropization time Tiso is likely located in some transition 
(isotropization) region. 

of the early-time behavior of the system to that of the CGC demonstrates that our conclusions 
may be relevant to heavy ion collisions in the real world. 

Fig. [2] summarizes our knowledge of the energy density in heavy ion collisions as a function 
of r for the strongly-coupled system considered above and in [1]. The energy density starts out 
as a constant, after which the system undergoes an isotropization transition (denoted by the 
dashed line in Fig. [2]) resulting in Bjorken hydrodynamics. The exact shape of the dashed curve 
is indeed unknown: what we plot in Fig. [2] is a simple interpolation between the early- and 
late-time regimes. More work is needed to develop a better understanding of this important 
intermediate region. 

Finally, above we estimated the time of the transition from our regime to that of Janik 
and Peschanski. The estimate is given in Eq. (IHTl) . If applied to RHIC, it appears to give the 
isotropization time Tiso ~ 0.3 fm/c, consistent with the results of hydrodynamic simulations 
[23-30]. 

In conclusion let us point out that the method of requiring the metric to be a real and 
single-valued function of coordinates could be successfully applied to other related problems. 
Recently o(l/r^/^) and o(l/r^/^) corrections to the solution of Eq. fl26l) have been calculated 
in [53, 54] for late proper times. One expands [53, 54] 

a{T,z) = a(w) + ai(w)— ^ 

and finds the subleading coefficients ai{v) and a2{v) by solving Einstein equations perturba- 
tively. [The same expansion is done for b{T, z) and c{t,z).] Taking the o(l/r^/^) correction to 
the coefficient a{v) in Eq. ( l26l) . denoted by 02 (f) and given in Eq. (20) of [53], one can see that 
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it also has a term proportional to ln(l — f ^/3), just like the solution in Eq. (126|1 with A = 



-4/3: 



1 3 71^ 

a2{v) = ( — -= ) ln(l — Dv^) + terms without branch cuts. (84) 



,4^3 

Here rjo is the proper-time independent coefficient related to shear viscosity r] by [53, 54] 

T] = — . (85) 
r 

Again we want the metric to be real and single-valued, hence the coefficient in front of ln(l — 
D v^) has to be an integer. Note however that a2{v) enters the metric with a prefactor of 1/r^/^ 
as can be seen from Eq. (!83|) . such that the coefficient in front of the logarithm in the metric 
varies with time. Therefore we conclude that the coefficient in front of ln(l — Dv^) can only 
be zero: this would insure that the metric is real at all times. Putting the coefficient in front 
of ln(l — Df^) in Eq. fl84p to zero yields 

V'o = ^ (86) 

in agreement with the result obtained in [53] using non-singularity of the curvature invariant 
criterion, and, as shown in [53], in agreement with the Kovtun-Policastro-Son-Starinets (KPSS) 
viscosity bound [43-47]! 

Here we have to caution the reader that similar branch cuts in other metric coefficients, 
like bi{v), b2{v) and Ci{v), C2{v) of [53] can not be eliminated this easily. (ai(w) has no branch 
cut singularity [53].) It is possible that resummation of the whole series of the type shown in 
Eq. (155]) for 6(r, z) and c(r, z) would be required to ehminate such branch cuts. We leave the 
investigation of this important question for future work. 
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